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Abstract 



We consider an anisotropic brane-world model with Bianchi type I and V geometry, without 
mirror symmetry or any form of junction conditions. The generaUzcd Chaplygin gas, which in- 
I terpolates between a high density relativistic era and a non-relativistic matter phase, is a popular 

candidate for the present accelerated expansion of the universe. Considering the generalized Chap- 
O ' lygin gas as a geometrical dark energy, we obtain the general solutions in an exact parametric form 

for both Bianchi type I and V space-times. Finally, we study the behavior of the observationally 
i important parameters such as the shear, anisotropic and deceleration parameter in this model. 



PACS numbers: 04.50.-h, 04.20.-q, 04.20.Jb 
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^ ■ 1 Introduction 
p 

' From a large number of observational evidence, the observable universe is presently undergoing an 

accelerated expansion. As an alternative to both the cosmological constant and quintessence, it is also 
] possible to explain the acceleration of the universe by introducing a cosmic fluid component with an 

exotic equation of state, called Chaplygin gas. The Chaplygin gas model describes a transition from 
a universe filled with dust-like matter to an accelerated expanding stage. The Generalized Chaplygin 
Gas (GCG) model, introduced in [1] and elaborated in [2], is described by a perfect fluid obeying an 



I exotic equation of state 



Pch = --IT' (1) 

Pch 



where ^ is a positive constant and < a < 1. The original Chaplygin gas corresponds to a = 1. 
However, recent analyses based on the latest of type-la Supernovae data have yielded rather surprising 
results, namely that a > 1 and that there is a degeneracy between the GCG and xCDM models in 
the form of a phantom-like energy component [30]. An attractive feature of the model is that it 
can naturally explain both dark energy and dark matter [3]. Within the framework of Friedmann- 
Robertson- Walker cosmology, this equation of state leads, after having inserted into the relativistic 
energy conservation equation, to an energy density evolving as 



Pch 



A + 



3{l+a) 



1 



l + Q, 

(2) 
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where a is the scale-factor of the universe and C is an integration constant which should be positive 
for a well-defined pch at all times. Hence, we see that pch ~ at early times, that is, pch behaves 
as matter while at late times it behaves like a cosmological constant pch ~ constant. The Chaplygin 
gas appears in the stabilization of branes in Schwarzschild anti-de Sitter (AdS) black hole bulks as a 
critical theory at the horizon [4] and in the string analysis of black holes in three dimensions [5] . The 
Chaplygin gas also appears as an effective fluid associated with d-branes [6] and can also be derived 
from the Born-Infeld action [7]. An interesting range of models have been found to be consistent 
with the SNe la data [8], CMB experiments [9] and other observational data [10]. The cosmological 
implications of the Chaplygin gas model have been intensively investigated in the literature [11, 12 , 13] . 

The idea that our familiar 4-dimensional (4D) space-time is a hypersurface (brane) in a 5D space- 
time (bulk) [14, 15, 16] has been under detailed elaboration during the last decade. According to 
this brane-world scenario, all matter and gauge interactions reside on the brane, while gravity can 
propagate in the 5D space-time. Several brane-world cosmologies have been proposed in the context 
of the Randall-Sundrum (RS) formulations [15], defined in a 5-dimensional anti-de Sitter space-time. 
The dynamics of these models feature boundary terms in the action and sometimes mirror symmetry, 
such that bulk gravitational waves interfere with the brane-world motion. This usually comes together 
with junction conditions producing an algebraic relationship between the extrinsic curvature and the 
confined matter [17, 18]. The consequence is that the Friedman equation acquires an additional term 
which is proportional to the square of energy density of the confined matter field [19, 20]. This term 
was initially considered as a possible solution to the accelerated expansion of the universe. However, 
soon it was realized to be incompatible with the big bang nucleosynthesis, requiring additional fixes 
[20]. 

Brane-world scenarios under more general conditions and still compatible with the brane-world 
program have also been rather extensively studied over the past decade where it has been shown that it 
is possible to find a richer set of cosmological solutions in accordance with the current observations [21] . 
Under these conditions, without using Z2 symmetry or without postulating any junction condition, 
Friedman equation is modified by a geometrical term which is defined in terms of the extrinsic 
curvature, leading to a geometrical interpretation for dark energy [22]. There has also been arguments 
concerning the uniqueness of the junction conditions. Indeed, other forms of junction conditions 
exist, so that different conditions may lead to different physical results [18]. Furthermore, these 
conditions cannot be used when more than one non-compact extra dimension is involved. Against 
this background, an interesting higher-dimensional model was introduced in [23] where particles are 
trapped on a 4D hypersurface by the action of a confining potential. The dynamics of test particles 
confined to a brane by the action of such potential at the classical and quantum levels were studied 
in [24]. In [25], the same brane-world model was studied, offering a geometrical explanation for the 
accelerated expansion of the universe. A geometrical explanation for the generalized Chaplygin gas 
was considered in [26] along the same line. In this paper, we consider an anisotropic brane-world 
model with Bianchi type I and V geometries filled with a geometrical Chaplygin gas. The behavior 
of the observationally important physical quantities is studied in this scenario. 

2 The model 

The embedding of the brane-world in the bulk plays an essential role in the covariant formulation of 
the brane-world gravity, because it tells us how the Einstein-Hilbert dynamics of the bulk is transferred 
to the brane. However, there are many different ways to embed a manifold into another, classified as 
local, global, isometric, conformal, rigid, deformable, analytic or differentiable. The choice of one or 
other depends on what the embedded manifold is supposed to do. 

Let us present a brief review of the model proposed in [22]. Consider the background manifold 
V4 isometrically embedded in a pseudo-Riemannian manifold Vm by the map y : V4 ^ Vm, with 
m = 4 + n components such that 

QAByiy!l=9^.u, QAByf^M^ = 0, GABM^Mi" = gab, (3) 
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where Qab {9^lv) is the metric of the bulk (brane) space Vm{V4) in arbitrary coordinates, {3^^} {{x^}) 
is the basis of the bulk (brane) and Af^ ^^rc n normal unit vectors orthogonal to the brane. According 
to Nash [27], we may continuously perturb along a normal direction in the bulk, to obtain another 
submanifold of the same bulk, provided the embedding functions remain regular. Perturbation of V4 
in a sufficiently small neighborhood of the brane along an arbitrary transverse direction ^ is given by 

2"^ = ^^ + {C^y)^, M^=M^ + (c^M)^ = AT^, (4) 

where £. represents the Lie derivative and (a = 1, 2, n) is a small parameter along M^, param- 
eterizing the extra noncompact dimensions. By choosing ^ orthogonal to the brane, we ensure gauge 
independency and have perturbations of the embedding along a single orthogonal extra direction Ma 
giving local coordinates of the perturbed brane as 

The above assumptions lead to the embedding equations of the perturbed geometry 

g^B^i^,^ = ffM-, QabZ^^M^ = g^a, QabM^M^ = gab. (6) 
From these equations it follows that 

g^^Zp;', = - g^'N^N^, (7) 
and also the components of the perturbed geometry 

Qixv = QabZ^Z^ = gn„ - 2^"'Knj,a + 9°''^ K^aaK^pb + g^'^^ij.ca^udb , (8) 

g^b = QabZ^^K = eA^ab, (9) 

gab = QAB^aJ^b = 9ab, (10) 



Kfxua = -QABM'a,ixZ^v = ^tiva - ^ 9 K^ccaKvfSb + 9"" ^ixca-^vdh , (H) 
A^ab = QabM^^^M^ = A^ab, (12) 

where A^ab represents the twisting vector fields and Kfj^i/a-, Kfj_ua represent the extrinsic curvature of 
the original and perturbed brane respectively. Comparing equations (8) and (11), we obtain 



\dg^^ 
2 ' 



K,ua = -^^4^, (13) 



which is the generalized York's relation and shows how the extrinsic curvatTirc propagates as a result 
of the propagation of the metric in the direction of extra dimensions. The components of the Riemann 
tensor of the bulk written in the embedding vielbein {Z^^N^}, lead to the Gauss-Codazzi and Ricci 
equations, respectively [28] 

Ral3-i5 = '^9"'^Ka[-yaKs]0b + T^ABCDZ^Z^Z^Z^g , (14) 

2-ft^Q[7c;6] = '^9°''' A[yacKs]ab + T^abcdZ^M^ Z^Z^ , (15) 

2^[7ab;5] = -'^9'''^ A^^^aAsjdb " 9^''^ K^^^aKsjdb " T^ABCoJ^a^I^ Z^Z^^ , (16) 
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where 'R-abcd and Ra/S'yS are the Riemann tensors for the bulk and the perturbed brane respectively. 
Contracting the Gauss equation (14) on a and 7 we find 

R^u = {Ki^acK - KcK^,'') + Uab^P^, - g'^^'TlABCDM^Zlz'^.M^. (17) 

A further contraction then gives the Ricci scalar 

R = {K^.aK^"'' - KaK'^) +n- 2g^''TlABN^K . (18) 
Therefore, the Einstein-Hilbert action for the bulk geometry in m-dimensions can be written as 

S = /"(7^-2A(''))V^c^'"x 

= ^/[R- {K^^aKi^-" - KaK'^) + Ig^^'TlABM^M^ - 2A('')] Vad'^x. (19) 

Variation of the action with respect to Qab gives the Einstein field equations in the bulk 

TIab - ^UGab = a*TXB - ^^^^Gab, (20) 

where a* = j^m--2 and A^''^ is the cosmological constant of the bulk space. The vielbein components 
of the energy-momentum tensor are given by 

T;, = TXsZp^„ t;, = t\bZ^^m^^ r:, = ti^aa^^a/^ . (21) 

The tangent components follow from the contractions of equation (20) with Z^Z^. After using 
equations (17) and (18) we obtain the "gravi-tensor" equation 

- \R9i.v - + g'"'TiABCDNtzlZ%N^ - g^^'TlABM^M^g^,, = a,T;, - A^'^g^,, (22) 

where 

Qm- = 5"' {KP^Kp,, - KaK^,,) 1 (K«/3a^"^" - KaK'^) g^^u. (23) 

By a direct calculation we can see that the extra term Q^y is an independently conserved quantity. 
On the other hand, the trace of the Codazzi equation gives the "gravi-vector" equation 

Finally, the "gravi-scalar" equation is obtained from equations (18) and (20) 



R - K^yaK'""'' + KaK'' = -2a. ( g^'x:, - J^T* ) + 2A(^) 



(n-l)(n + 4) 

n — 



(n + 2) 



(25) 



In its most general form, without assuming extra dimensional matter, the confinement hypothesis 
states that the only non-vanishing components of Tab are the tangent components T^i, representing 
the confined sources. Therefore we set 

Oi*T*y = 8TrGTfj_„, a*T*a = 0, a^T*i, = 0. (26) 

Equations (22)-(25) represent the most general equations of motion of a brane-world, compatible with 
the differentiable embedding in a m-dimensional bulk defined by the Einstein equations. Clearly, the 
usual Einstein equations are recovered when all elements of the extrinsic geometry are removed from 
those equations. 

The geometrical approach considered here is based on three basic postulates, namely, the confine- 
ment of the standard gauge interactions to the brane, the existence of quantum gravity in the bulk 
and finally, the embedding of the brane-world. All other model dependent properties such as warped 
metric, mirror symmetries, radion or extra scalar fields, fine tuning parameters like the tension of the 
brane and the choice of a junction condition are left out as much as possible in our calculations [21]. 
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3 Field equations and observational parameters in anisotropic brane 



In what follows, we will investigate the infincncc of the extrinsic curvature terms on the anisotropic 
universe described by Bianchi type I and V geometries. From a formal point of view these two 
geometries are described by the line element 

ds'^ = -dt^ + al{t)dx'^ + ai(i)e-2^^dy2 ^ ^2^^^y-2px ^^2 ^ ^27) 

where ai{t),i = 1,2,3 arc the expansion factors in different spatial directions. The metric for the 
Bianchi type I geometry corresponds to the case /3 = 0, while for the Bianchi type V we have /3 = 1. 
Let us define the following variables [29] 

3 . 3 

v = l[ai, Hi = ^,i = 1,2,3, 3H = Y,Hi, AHi = Hi - H,i = 1,2,3. (28) 

i=i i=i 

In above equations, v is the volume scale factor, Hi,i = 1,2,3 arc the directional Hubble parame- 
ters, and H is the mean Hubble parameter. The physical quantities of observational importance in 
cosmology are the expansion scalar 0, the mean anisotropy parameter A, the shear scalar parameter 
a^, and the deceleration parameter q, which are defined according to 

e^sH, 3^=^:(^)^ .^4.,,.«=l^:i^?-3i^^ ,=1(1) -1. (29) 

1=1 1=1 

The sign of the deceleration parameter indicates how the universe expands. A positive sign for q 
corresponds to the standard decelerating models whereas a negative sign indicates an accelerating 
expansion at late times. We note that A = for an isotropic expansion. 

Let us assume that the confined source on the brane is a perfect fluid with a linear barotropic 
equation of state, namely p = — l)p with 1 < 7 < 2. In this paper we restrict our analysis to a 
five-dimensional bulk with a constant curvature characterized by the Riemann tensor 

T^ABCD = hiGAcGsD — GadGbc), (30) 

where A;* denotes the bulk constant curvature. In the flat case A;* = and in the de Sitter and anti-de 
Sitter cases we may write A;* = respectively. Assuming g^^ = 1 and using the Gauss-Codazzi 

equations, we obtain 

Rap-yS = {Ka-yKps — KasKp^) + K{ga-ygp5 " gaSgp-y), (31) 



The equations of motion derived in the previous section can be obtained directly from equations (31) 
and (32). The result is Einstein equations as modified by the presence of the extrinsic curvature 

Gfj,i, = SirGTfj^i, - kg^^ + Q^^, (33) 

where A = — 3/c* -|- A*^^^ is the effective cosmological constant in four dimensions with Q/ti^ being a 
completely geometrical quantity given by 

Q^, = {KK^, - K^aK) + 1 (k^^K'^'^ - K^) g^„ (34) 
where K = g'^'^K^^. Using the York relation 

K^iua = —^-^, (35) 
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we realize that in a diagonal metric, K^j^a is diagonal. After separating the spatial components, the 
Codazzi equations reduce to (here a, 7, cr = 1, 2, 3) 

i^V,, + K%j:-p, = K\,^^ + K/^,„rv, (36) 



i^V,o + ^K%a = ^.5%K\,, z = 1, 2, 3. (37) 



The first equation gives K^ia,a = for cr 7^ 1, since K^i^ ^^^^ not depend on the spatial coordinates. 

Repeating the same procedure for 0,7 = = 2, 3, we obtain K'^2a,a — a ^ 2 and ^^3^ ^ = 
for (7 7^ 3. This shows that K^ia-> functions of t and the choice K^^^ = = 

if^3„ = ba{t), where ba{t) are arbitrary functions of t, would simplify our analysis. Now from the 
second equation we obtain 

4 + -ba = -i^V> i = l,2,3. (38) 
Summing equations (38) we find 

KoOa = -(^^+bay (39) 

For 11,1' = 1, 2, 3 we obtain 

K^va = baQnu- (40) 

Note that in equation (33) we have considered a five-dimensional bulk (m = 5) , thus the functions 
(a = 1,2, ...,m — 4) reduce to only one function, namely 61. Denoting bi = b, 6 = ^ and = ^, we 
find from equation (34) that 

(29 \ 
-Q+lj^/.., /i,i/=l,2,3, goo = 362. (41) 

As we have noted before, Q^i, is an independently conserved quantity, that is Q^^^ = 0, suggesting 
an analogy with the energy momentum of an uncoupled non-conventional energy source. We see 
that solution (41) depends on the arbitrary function b{t). To find the dynamical role of this function 
and to compare the compatibility of such geometrical model with the present experimental data, we 
assume Q^i, to be a conserved energy-momentum tensor and take the GCG as an example 

1 A 

QlJ-i' ^ 3—77 [{Pch + Pch)UiiUu + Pch9iJ.iy] , Pch = a' (^2) 

87rGr p^f^ 

where A and a arc positive constants. Comparing Qfj,u, fJ-, i' = 1,2,3 and Qoo from equation (42) with 
the components of and Qqq given by equation (41), we obtain 

362 \ 362 

Use of the above equations leads to an equation for b{t) 



for which the solution is 



/ 8ttg \ 2 r , . c ^ 2(i+a) 
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where C is an integration constant. Using equation (43) and this solution, the energy density of the 
GCG becomes 



Pch = Pcho 



(l+a) 



(46) 



where pch^ is the GCG density at the present time, A, 
to the speed of sound for the GCG today, v1 = aAp^j^"^^"^^ and C 



^/'c/io^^"^ is a dimensionless quantity related 



(l+o) 
Pcho 



A. From equation (46), 



we see that for = the GCG behaves like matter, whereas for = 1 it behaves as a cosmological 
constant and when < As < 1, the model predicts a behavior as that between a matter phase in the 
past and a negative dark energy regime at late times. This particular behavior of the GCG inspired 
some authors to propose a unified scheme for the cosmological "dark sector," an interesting idea 
which has been considered in many different contexts. 

Now, using the geometrical energy density for Q^^, the field equations on the anisotropic brane 
become 



3H + Y,Hf = A- 47rGpo(37 - 2)^;"^ + SirGpcho 



As + 



yil+a) 



As 



(47) 



~ivH) = 2/^2^-2/3 + A - 47rGpo(7 - ^)v-'' + SnGpcho 



As 



(1 - As, 

,,(l+a) 



l+a 



As + 



2^(1+") 



.(48) 



For /? = we obtain the field equations for Bianchi type I geometry, while (3=1 gives the Bianchi 
type V equations on the anisotropic brane. Using the relation = we can rewrite equation (48) 
in the form 



V = 6/3\^/^ + 3Av - 127rGpo(7 - 2)?^^"^ + 247TGvpcho 
The general solution of equation (49) becomes 
t-to-- 



As + 



{I -As 



l + a 



As + 



:i - A, 



2i;(i+") 



dv 



^9/?2t;4/3 + 3A^2 + 24TrGpov^-^ + f{v) + c' 



where 



f{v) 



247rGpcho 
2 + a 



V\\ + 



Asv 



1+q\ l+o 



I -As 



1 + As (v^+" - l) 



l+a 



\l + a l + a l + a Ag — 1 



+ (2 + a)(l-^)2Fi 



a 2 + a AsV^+'' 



1 + a'l + al + a' Ag 



(49) 



(50) 



(51) 



and C is a constant of integration. The time variation of the physically important parameters de- 
scribed above in the exact parametric form, with v taken as a parameter, is given by 



e = 3H 



9^2^4/3 + 3A^;2 + 247rGpov^-^ + f{v) + C 



(52) 



o-i = aoiV^^^ exp 



/ 



hidv 



v^QP^v"^/^ + 3Av^ + 247rGpov2-7 + f{v) + C 



z = l,2,3. 



(53) 
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9/52^,4/3 + 3At,2 + 247rG/9oi;2-7 + /(i;) + C' ^ ' 



„{l+a) 



1+a 



-I- 2t,(l+a) 



3^2^4/3 + Ai;2 + STrGpot^^-^ + ^/(t^) + C 



(56) 



where /ij, i = 1, 2, 3 are constants of integration and h"^ = Yl\=i The GCG provides an interesting 
candidate for the present accelerated expansion of the universe without resorting to an effective 
cosmological constant. To this end we consider A = and show that, within the context of the 
present model, the geometrical Chaplygin gas can be used to account for the accelerated expansion 
of the universe. Therefore for A = and a = 1, the dynamics of the Bianchi type I universe is 
controlled by the Chaplygin gas equation of state parameter As. In this case, the expansion, scalar 
factor, anisotropy, shear and deceleration parameters are respectively given by 



Q = 2,H 



247rGpot;2-7 + 247rGpe^„^;2[^, + (1 - ^)t;-2]V2 + c 



(57) 



aoiV^^^ exp 



hjdv 



^247rGpoi;4-7 + IA-kGp^HovMs + (1 - ^)^;-2]V2 + Cv'^ _ 



1,2,3, (58) 



^_ 3^i2 

" 247rGpoi;2-7 + 247rGpeho^^'[^ + (1 - As)v-^Y'^ + C 



(59) 



(60) 



127rGpo(2 - 7)z;2-7 + 2A7:GpchoV 



-1/2 



87rGpoi;2-7 + 87rGpch,y[As + (1 - As)v-^]y^ + C 



(61) 



In figure 1 we present the dynamics of the deceleration parameter for different values of Ag and for 
7 = 1 and a = 1. This behavior is much dependent on the range of the values that As can take. For 
having an accelerating universe, the value of Ag should lie in the range < < 1. In the initial 
stage the evolution of the Bianchi type I brane universe is non-inflationary, but in the late time limit 
the brane universe ends in an accelerating stage. 

For a better understanding of the behavior of the mean anisotropy parameter, let us consider it 
as a function of the volume scale factor 



A{v) 



3^2 



9^2^4/3 + 3Ay2 + 247rGpo?;2-7 + f{v) + C" 



(62) 



where f{v) is defined by equation (51). The behavior of the anisotropy parameter at the initial state 
depends on the values of a and Ag. For original Chaplygin gas with a = 1, we obtain an accelerating 
universe when < < 1. From equation (62), in the limit v ^ and taking < < 1, we find 
that the initial state is always anisotropic A{v) ^ 0. The behavior of the mean anisotropy parameter 
of the Bianchi type I and V geometries is illustrated, for a = 1 and different values of Ag, in figure 2. 
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Figure 1: Left, deceleration parameter q for the geometrical Chaplygin gas filled Bianchi type I brano universe as a 
function of time and right, the same parameter in the Bianchi type V brane universe for 7 = 4/3, a = 1 and As = 0.3 
(solid line), = 0.5 (dashed line). As = 0.8 (dot-dashed line) with A = 0. 



A(t) A(t) 




Figure 2: Loft, anisotropy parameter A for the geometrical Chaplygin gas filled Bianchi type I brane universe as a 
function of time and right, the same parameter in the Bianchi type V brane universe for 7 = 4/3, a = 1 and As = 0.3 
(solid line). As = 0.5 (dashed line). As = 0.8 (dot-dashed line) with A = 0. 




Figure 3: Left, shear scalar ct^ for the geometrical Chaplygin gas filled Bianchi type I brane universe as a function of 
time and right, the same parameter in the Bianchi type V brane universe for 7 = 4/3, a = 1 and As = 0.3 (solid line), 
As = 0.5 (dashed line). As = 0.8 (dot-dashed line) with A = 0. 
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Figure 4: Left, deceleration parameter q of the Bianchi type I brane universe as a function of time and right, anisotropy 
parameter A of the Bianchi type I brane universe as a function of time for 7 = 4/3, a = 3.75 and = 0.3 (sohd line), 
As = 0.6 (dashed line), A^ = 0.936 (dot-dashed hne) with A = 0. 

The behavior of this parameter shows that the universe starts from a singular state with maximum 
anisotropy and ends up in an isotropic de Sitter inflationary phase at late times. The time variation 
of the shear parameter is represented for different values of Ag, in figure 3. 

It is worth mentioning that although most of the works on GCG cosmology have assumed a value 
for a compatible with < a < 1, it has been shown that the type-la Supernovae data favors a > 1 
[30, 31, 32]. The best fitted values suggested for a > 1 are a = 3.75 and As = 0.936 [30]. In figure 
4 we have plotted the deceleration and anisotropy parameters of the Bianchi type I geometry for 
a = 3.75 and different values of Ag. This behavior shows that for a > 1 with < < 1, the 
geometrical model presented in this work is in agreement with observational data. 

4 Conclusions 

In this paper we have shown that dark energy may be considered as a consequence of the extrinsic 
curvature in a brane-world scenario and extended the predictions of the geometrical matter in the 
more general case when the relation between p and p is not linear [26] . In this work, we have studied 
the Bianchi type I and V geometries, seen as a brane-world embedded in a five-dimensional bulk of 
constant curvature, without Z2 symmetry or any form of junction condition. We have shown that the 
geometrical Chaplygin gas may be used to account for the accelerated expansion of an anisotropic 
universe. We have also obtained the general solutions in an exact parametric form for both Bianchi 
type I and V geometries and studied the behavior of the observationally important parameters. 

The study of anisotropic homogeneous brane-world cosmological models has shown that in the 
framework of the RS models, brane universes are born into an isotropic state [11], whereas in the 
model presented here the universe starts as a singular state with maximum anisotropy and reaches 
an isotropic state in the late time limit, in agreement with the standard AD cosmology. 
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